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ABSTRACT 


The results of an investigation of the influence of simulated 
turbulence on aircraft handling qualities is presented. Pilot opinions 
of the handling qualities of a light general aviation aircraft were 
evaluated in a motion— base simulator using a simulated turbulence 
environment. A realistic representation of turbulence disturbances is 
described in terms of rms intensity and scale length and their random 
variations with time. The time histories generated by the proposed 
turbulence models showed characteristics which appear to be more similar 
to real turbulence than the frequently-used Gaussian turbulence model. 

In addition, the proposed turbulence models can flexibly accommodate 
changes in atmospheric conditions and be easily implemented in flight 
simulator studies. 

Six turbulence time histories, including the conventional Gaussian 
model, were used in an IFR-tracking and a landing approach task. The 
realism of each of the turbulence models and the handling qualities of 
the simulated airplane were evaluated. Analyses of pilot opinions show 
that at approximately the same rms intensities of turbulence, the handling 
quality ratings transit from the satisfactory level, for the simple 
Gaussian model, to an unacceptable level for more realistic and composite- 
ly structured turbulence models. 
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SECTION I 


INTRODUCTION 

Simulated time histories of aircraft motion in a turbulence 
environment are required in a variety of engineering applications, 
and their use appears to be increasing as more intricate and 
sophisticated design studies are attempted. As an example, the use 
of flight simulators for the study of airplane handling and ride 
quality has proven to be more valuable when disturbances in the form 
of artificially simulated turbulence are introduced into the system. 
Several methods have been used to generate turbulence signals; each 
one aimed at realizing the actual atmosphere as closely as possible, 

A realistic representation of turbulence becomes especially important 
in the simulation of future aircraft with high sensitivity to 
turbulence, as even light to moderate turbulence may seriously degrade 
their controllability and ride quality. Low altitude atmospheric 
turbulence critically affects the evaluation of vehicle handling 
qualities, pilot work load, ride quality, and other design factors. 
Several empirical studies Q->2,3) have shown that low altitude, clear 
air atmospheric turbulence is only locally isotropic, i.e, , isotropic 
over a finite range of scale lengths. The proposed gust model accounts 
for the anisotropy of typical low altitude clear air turbulence by 
randomly varying the rms velocities and scale length of the gust field. 
The scale lengths predicted by either the Von Karman or the Dryden 
models C4-) are large compared to real atmospheric turbulence and hence 
the scale length distribution is modified to achieve compatibility. 
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With a suitable combination of scale length and intensity 
distribution, the proposed model will simulate various atmospheric 
conditions characterized by altitude, stability, and terrain. This 
new model is mechanized to be included in a flight simulator experiment 
in order to determine to what extent the pilots are sensitive to 
changes in atmospheric conditions and the realism of the model. The 
flight simulator experiment is conducted for two sets of landing 
conditions. The first requires the pilot to follow an IFR-tracking 
task with no out-the-window cues provided. In the second condition 
the simulator is equipped with a visual display which provides a 
realistic landing approach scene. The following chapters describe 
the proposed turbulence model and the flight simulator experiment 
in detail. 




SECTION II 


LITERATURE SURVEY 

In this chapter statistical properties of atmospheric turbulence 
are reviewed and presently-used simulation techniques are discussed. 

A review of basic definitions in probability and statistics is 
included in Appendix A. 

2.1 Properties of Atmospheric Turbulence 

Simulations of aircraft flying through atmospheric turbulence 
require a realistic model of the physical environment. Therefore, 
simulation studies in general begin with a study of the real 
atmosphere. In references 1-3, atmospheric data have been reported 
characterizing various atmospheric conditions for variation in 
terrain, stability, altitude, temperature, time, season, and 
geographic location. This data has been suitably modified to 
establish a basis of comparison for the simulated turbulence field. 

The following criteria are used as the bases of comparison: 

a) Output Statistics 

Mean and standard deviations of the gust velocities. 

b) Probability Distribution 
Cumulative probability 
Probability density 

Fourth and sixth normalized moments 

c) Patchiness of the Field 

d) Power Spectral Density 

e) Element of Surprise 
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Each of these properties will be discussed from the standpoint of real 
atmospheric turbulence. 

Mean ; Analysis of several sets of data presented in Reference (1) 
indicates that the mean velocity of atmospheric turbulence is 0.0 + 0*1 
ft/sec CO +_ 0.03 m/sec) . 

Standard Deviation ; The standard deviation of the velocity field 
for low altitude clear air turbulence is 3.0+1.31 ft/sec CO. 91 + 0.4 
m/sec) , Typical values for various conditions are listed in Reference 
C4) as ; 

0^-2 ft/sec CO. 61 m/sec) for light turbulence 

a=4 ft/sec (1.22 m/sec) for moderate turbulence 
u 

o^ = 6 ft/sec (1.82 m/sec) for severe turbulence 
where u is the longitudinal gust component. 

Probability Distribution ; The probability distribution of a 
random process provides information concerning the range of values 
assumed by that function and the frequency with which they occur. 

As there is little experimental data available which distinguishes 
between probability distributions of different gust components, no 
distinction will be made here. 

Probability Density Distribution ; Figure 1 presents data from 
Reference (5) showing a typical probability density distribution of 
atmospheric turbulence velocity. The departure from the Gaussian 
curve clearly indicates increased probabilities of large and small 
gus ts , 

Normalized Central Moments ; The fourth and sixth moments of low 
altitude real atmospheric turbulence are M^ = 3,5 and M^ 
respectively. (5) 


= 21.7, 




Probability % 
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w - ft/sec (m/sec) 


Figure 1 Typical Probability Density Function of Atmospheric 
Turbulence 
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Patchiness : It is known that turbulence has a non-Gaussian 

patchy structure which seems to occur in bursts of relatively intense 
motion separated by areas of relative calm. Figure 2 shows typical 
patchy characteristics for a 40-sec sample of real atmospheric 
turbulence , 

Power Spectral Density (PSD) : The PSD of a random process 

provides information on the average contribution to the process from 

the frequency components which make it up. Figure 3 presents a 

typical plot of PSD of low altitude clear air atmospheric turbulence. 

It may be observed that at high frequencies, the spectral density 

-2 

varies as inverse square of frequency (w ) . On the other hand, at 
low frequencies the PSD is characterized by a horizontal asymptote. 
Two convenient mathematical forms are used to represent the power 
spectra of atmospheric turbulence. These are: 
a) Von Karman Spectra 
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b) Dryden Spectra 
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Ug = initial total velocity 
till 

L = scale for i turbulence velocity 

i 

cr^ " rms gust intensities 
w = frequency 

i = u, v, w gust components. 

The Von Karmen spectral shapes, although accurate, are not con- 
venient for turbulence modelling work since they cannot be matched 
using linear filters. This is due to the noninteger power appearing 
in the denominators. Thus, in order to avoid computational complexity 
in this report, the Dryden form is adopted. 

Element of Surprise : More often than not, real atmospheric 

turbulence, when encountered, presents an element of surprise. It is 
not easy to formulate a model of this phenomenon in terms applicable 
to flight simulator work. It seems that a measurement of ’’sudden jump” 
in the velocity field can be used as a possible criterion to describe 
this phenomenon. Relative frequency of "sudden jump” of atmospheric 
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turbulence can be compared to the simulated turbulence field. Changes in 
aircraft orientation angles can also be used to measure this phenomenon. 
2,2 Presently-Used Simulation Techniques 

In this section several pvesently-used simulation techniques are 
discussed from the standpoint of their statistical realism and suit- 
ability for use in flight simulators . 

Measured Turbulence Field ; Flight recordings of atmospheric 
turbulence is perhaps the most obvious method of producing a realistic 
simulation. There can he little argument as to whether or not these 
time histories are an accurate and realistic representation. However, 
it is difficult to adjust the measured time histories to allow for 
conditions other than those for which it was recorded. No allowances 
can be made for changes of altitude or different atmospheric conditions. 
Another serious drawback is that the recorded time histories are fixed 
in length. Extended run times, therefore, cannot be accommodated 
without repetition. From the simulation point of view, the pilots 
tend to recognize some of the characteristics of the turbulence field 
and develop an intuition for predicting the field. This defeats the 
purpose of an artificially simulated turbulence field, which is to 
provide unpredictable external disturbances. It can, therefore, be 
concluded that flight recordings of atmospheric turbulence are not 
suitable for the simulation of typical turbulence. 

Sum of Sine Waves : Reference C5) describes this method in summary 

form. This technique involves superimposing several sinusoidal waves 
of different frequencies and amplitudes. The resultant is used to 
represent time histories of turbulence. One obvious disadvantage of 
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this method is that it contains only a finite range of frequencies 
whereas actual atmospheric turbulence consists of an infinite number 
of frequency components. 

Results of this simulation are not available but the model can 
justifiably be discarded on the basis of its inadequacy in matching 
the frequency content. 

Method of Orthogonal Functions ; In this method (7) , the recorded 
time histories of turbulence are decomposed into eigenfunctions of a 
covariance matrix. The probabilistic structure of the eigenfunction, 
and the coefficients of each of the time histories are studied. 

Simulated time histories are then regenerated by suitably modifying 
the distribution of the coefficients. The available preliminary 
results show that this technique adequately models the frequency contents 
and also presents an element of surprise. However, this model fails to 
show a patchy non-Gaussian characteristic which is typical of the real 
atmosphere. In addition to the mathematical complexity of the 
technique, its application is limited since recorded time histories 
are needed. 

Gaussian Turbulence Model ; The classical method, most widely used 
for turbulence simulation, is the linearly filtered white noise tech- 
nique. Here the turbulence gust field is produced by passing white 
noise through a linear filter as shown in Figure 4a. The resultant 
signal is shaped so that the power spectrum and rms intensities match 
those of real turbulence, A Dry^en or Von Karman form (6) are normally 
used to model the power spectrum. This model is remarkably easy to 
implement and can be adjusted for any general power spectrum. However, 
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Component 


Figure 4a Gaussian Turbulence Simulation 



Figure 4b Derivative of Vertical Velocity 
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this model too falls short of reproducing the non-Gaussian patchy 
nature of real turbulence. Figure 4b compares the artificially 
simulated gust field using the Gaussian model (with a Dryden spectrum) 
and real atmospheric turbulence. It may be observed that the intensity 
for the Gaussian model is nearly constant whereas measured ("real") 
turbulence exhibits a patchy nature or intensity bursts. Test pilots, 
when exposed to this model in a flight simulator, rated the realism 
fair to poor. (5,6) 

Kon-Gaussian Turbulence Model ; Reference (6) presents a non- 
Gaussian turbulence model. Time histories are generated by multiplying 
too independent random variables, one to represent the turbulence 
within a patch and the other to represent the variation of intensity 
with time. Figure 5a shows two independent Gaussian white noise 
generators and linear filters, which produce Gaussian random variables, 
a(t) and b (t) . These variables are then multiplied to produce gust 
time histories. 

The non-Gaussian model proposed in Reference (5), a modification 
of the above, is shown in Figure 5b, Here a(t), b(t), and d(t) are 
independent Gaussian processes. The process c(t) is generated by 
multiplying a(t) and b(t). The resultant process, c(t), a modified 
Bessel process, is summed with d(t) to form the output, u(t) , The 
most remarkable achievement of this model is that the patchy character- 
istic and several statistical parameters of the simulated turbulence 
field can be varied simultaneously by varying the standard deviation 
ratio (R = cr /a,). However, when R is varied to achieve one set of 

C Q 

statistical properties, several other statistical parameters of 
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Figure 5a Non-Gaussian Turbulence Model (Reference 6) 
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Figure 5b Non-Gaussian Turbulence Simulation {Reference 5) 









16 

interest do not match real turbulence. In addition, due to the 
mathematical complexity, the mechanisation of this model on a flight 
simulator is complicated and expensive. 

It can be observed from the review of presently-used simulation 
techniques that there is a need for a new model which adequately 
matches real atmospheric turbulence and is simple to implement in 
flight simulator studies. None of the preceding models have the 
flexibility of simulating various atmospheric conditions characterized 
by altitude, stability, and terrain. It is, therefore, necessary to 
introduce a new turbulence model which is realistic and can flexibly 
accommodate changes in atmospheric conditions and be easily implemented 
in flight simulator studies. 



SECTION III 


PROPOSED GUST MODELS 


Of the simulation techniques described, the Gaussian turbulence 
model is the simplest to implement and least expensive computationally. 
The proposed turbulence models, modification of the Gaussian simulation 
technique, retain the simplicity of the Gaussian technique while 
adequately modelling the characteristics of real atmospheric 
turbulence. In this report three basic models are proposed: 

1) Modified Gaussian Model 

2) Rayleigh Model 

3) Variable Length and Intensity (VLI) Turbulence Model. 

3.1 Modified Gaussian Model 

A. block diagram of the modified Gaussian model is presented in 
Figure 6, Gaussian white noise, is passed through a linear filter, 
G^Cs) i = u, v, w, whose power spectrum is given by a Dryden model 
(e.g,, Eqs. 2,1.4 to 2.1,6). The mathematical form of linear filter 
G.(s) is given as follows: 



G 

v 


GO 



(3.1.1) 


(3.1.2) 
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(3.1.3) 


where ^ is the white noise power spectrum. 

The linear filter , described above, is modified to include random 
variation of rms intensities. Random numbers generated by A are 
passed through a distribution modifier to generate rms intensities. 

Time histories are then generated by passing Gaussian white noise, <{)q, 
through the linear filter modified by the distribution modifier. 

The patchy nature of atmospheric turbulence suggests that the 
turbulence field is composed of two components. One to represent 
variation of intensity within a patch and the other to represent 
variation of intensity with time (or from patch to patch) . The 
distribution modifier in this model, essentially, represents the 
variation of intensity with time. The level of turbulence within 
each patch is controlled by the magnitude of the rms intensity. 

The distribution modifier is the probability density function of 
the rms intensity. Analysis of several sets of atmospheric data 
characterized by various atmospheric conditions show that a truncated 
Gaussian distribution best fits the probability density of rms 
intensity. (1) 

rms Distribution Modifier; 

< m 2*1 

(3.1.4) 
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P - probability density function 
= rms intensity 

S = root mean square of rms intensity 
m ~ mean of rms intensity 
i = u>, v, w gust components. 

Equation 3.1*4 is completely described by the mean, m, and the root 
mean square, S, of the rms intensity, These variables have been 
derived from the data presented in Reference (1) characterized by 

(j 

terrain, altitude, and atmospheric stability, Table 1 represents the 
distribution modifier for two sets of atmospheric conditions. Through- 
out this report, the turbulence generated by these too distribution 
modifiers will be referred to as Model 2 and 3 (.Model 1 is Gaussian 
turbulence simulation) » 

The scale lengths for these models are given by the Dryden form: 

= h for h ^ 1750 ft (533,4 m) (3.1.5) 

L =1 =145 h 1/3 for h < 1750 ft (533.4 m) (3.1.6) 

u v 

L - h (3.1.7) 

w 

where h is the altitude. 

3,2 Rayleigh Model 

The Rayleigh model is derived from the modified Gaussian model by 
replacing the distribution modifier by a Rayleigh probability density 
function. The Rayleigh probability density function for rms vertical 
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TABLE 1 

DISTRIBUTION MODIFIERS (rms INTENSITY) 


Mean Variance 


rms Distribution Modifier 
Model 2 


Altitude: 

250 ft C76.2 m) 

a ft/sec 
u 

Cm/sec) 

3.1 

CO. 94) 

1.2 

CO. 37) 

Atmospheric 

Stability: 

Tins table 

a ft/sec 

V 

Cm/ sec) 

3.2 

CO. 97) 

1.2 

CO. 37) 

Terrain: 

Plains 

cr ft/sec 
w 

Cm/ sec) 

2.8 

CO .85) 

0.9 

(0.27) 


rms Distribution Modifier 
Model 3 


Altitude: 

750 ft (228 ,j 6 ra) 

a 

u 

ft/sec 

Cm/ sec) 

3.2 

CO. 97) 

0.8 

(0 . 24) 

Atmospheric 

Stability.: 

Unstable 

0 

V 

ft/sec 

Cm/ sec) 

3.5 

Cl. 07) 

1.0 

CO. 30) 

Terrain: 

Mountain 

a 

w 

ft/sec 

Cm/ sec) 

4.1 

Cl, 25) 

0.9 

CO. 27) 
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turbulence intensity is, is given by 


?Ca ) 

w 



exp 


c- 


1 

2 



C3.2.1) 


where C 2 is one-half the expected value of cr 2 , 

w 

Using Dryden spectrum models of real atmospheric turbulence, the 
value of C has been estimated in Reference 0*) to be 2.3 ft/sec (0-70 
m/iiec) . 

The rms intensity of the longitudinal, u, and the lateral, v, gust 
components are obtained from the relation: 



C3.2.2) 


The scale lengths are given by Equations 3.1.5 to 3,1.7. This will be 
referred to as Model 4. 

3.3 Variable Length and Intensity (VLI) Turbulence Model 

The VLI turbulence model includes, in addition to the rms 
distribution modifier, a scale length modifier. A block diagram of 
this model is presented in Figure 7. In addition to controlling the 
patchiness of the turbulence field, the time variations of scale length 
achieves numerical compatibility with the real atmosphere and further 
randomises the simulation. 

The scale length distribution modifier is derived from data 
collected in the LO-LQ-GAT Program Cl) for various combinations of 
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altitude, terrain, and atmospheric stability, Figures 8 and 9 show 
the fitted Gaussian distribution of scale length modifier for two sets 
of atmospheric conditions. The scale length distribution modifier is 
■assumed to have the form 


KL.) 


1 


exp 

s/2tF 



where 

P = probability density function 
= scale length of i^ component 
S = root mean square of scale length 
m ~ mean of scale length distribution 
i = u, v, w gust components 


(3,3.1) 


Table 2 presents the root mean square and mean of scale length 
distribution along with the rms distribution modifier for specific 
atmospheric conditions. The turbulence signal generated by these two 
atmospheric conditions will be referred to as Models 5 and 6. 



Probability (1CT 3 ) 


p Longitudinal Component 
O Lateral Component 

X Vertical Component 


Model : 5 

Altitude: 250 
Terrain: Plains 
Atm. Stability: Unstable 



Symbols represent empirical data.- 


Figure 8 Scale Length Distribution 
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•f- 

□ Longitudinal Component 
O Lateral Component 



4 . 

Symbols represent empirical data. 


Figure 9 Scale Length Distribution 






TABLE 2 


DISTRIBUTION MODIFIERS (SCALE LENGTH) 


Mean Variance 


ras Distribution Modifier 
Model 5 








Altitude: 250 ft {76,2 m) 

a 

u 

ft/sec 

(m/sec) 

3.1 

(0.94) 

1,2 

(0,37) 

Atmospheric 

Stability: Unstable 

a 

V 

ft/ sec 

(m/sec) 

3,2 

(0,97) 

1,2 

(0,37) 

Terrain: Plains 

cr 

w 

ft/sec 

(m/sec) 

2,8 

(0.85) 

0,9 

(0,27) 


L 

u 

ft (m) 


415 

(126,4) 

110 

(33,5) 

Scale Length Modifier 
Model 5 

L 

V 

ft (m) 


325 

(99.1) 

86,6 

(26.4) 


L 

w 

ft (m) 


335 

(102.1) 

83,1 

(25,3) 

rms Distribution Modifier 
Model 6 








Altitude: 750 ft (228,6 m) 

a 

u 

f t/sec 

(m/sec) 

3.2 

(0.97) 

0,8 

(0.24) 

Atmospheric 

Stability: Unstable 

a 

V 

f t/sec 

(m/sec) 

3,5 

(1.07) 

1.0 

(0.30) 

Terrain: Mountains 

a 

w 

ft/sec 

(m/sec) 

4,1 

(1.25) 

0,9 

(0,27) 


L 

u 

ft (m) 


^ ; 415 

(126,4) 

116.5 

(35.5) 

Scale Length Modifier 
Model 6 

L 

V 

ft (m) 


460 

(140,2) 

126,5 

(38,6) 


L 

w 

ft (m) 


425 

(129.5) 

132,9 

(40.5) 


SECTION IV 


THEORETICAL ANALYSIS OF PROPOSED MODELS 

In this section results obtained by statistical analysis of the 
gust velocity components for each of the six models will be discussed 
and compared with the properties of real atmospheric turbulence where 
possible. The statistical results have been obtained in the form of: 

1) mean and standard deviations 

2) normalized fourth and sixth moments 

3) probability density functions 

4) power spectral densities 

5) patchiness 

6) frequency of element of surprise. 

Table 3 tabulates the mean and standard deviation of gust 
components for each of the six models. It may be observed that the 
standard deviation varies from 2.6 to 5.2 ft/sec (0.79 to 1,58 m/sec) 
which is typical of low altitude clear air turbulence. 

Fourth and sixth moment characteristics are tabulated in Table 4. 
Within the limits of experimental error these characteristics for the 
VLI models are in fairly good agreement with the real atmospheric data 
obtained in Reference (5) . 

Since the cumulative probability and the probability density 
function essentially contains identical information, only the 
probability density function will be analyzed. Figures 10 to 15 are 
plots of probability density functions for the simulated cases. In 
order to compare these with real atmospheric turbulence, a Gaussian 
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TABLE 3 


MEAN AND STANDARD DEVIATION OF GUST COMPONENTS 
CIO Tain, sample) 


Gust Component 


Model 

Output 

u ft/sec 

v ft/sec 

w ft/sec 

Simulation 

No. 

Statistics 

(m/sec) 

Cm/sec) 

Cm/ sec) 

Technique 

1 

Mean 

0,08 

CO. 02) 

0.06 

CO, 02) 

-0.03 

C-0.009) 

Gaussian 

_L 

St. Deviation 

3.97 

Cl. 21) 

3.90 

Cl. 18) 

4,43 

Cl. 35) 

n 

Mean 

0,83 

CO. 25) 

-0,32 

C-0,09) 

-0.15 

C-0.04) 

Modified 

Z 

St, Deviation 

3.90 

Cl. 18) 

3,50 

Cl. 06) 

2.60 

CO. 79 

Gaussian 


Mean 

0,88 

CO. 27) 

-0.40 

C-0.12) 

0.06 

CO. 02) 

Modified 

J 

St. Deviation 

3,90 

Cl. 18) 

3.90 

Cl. 18) 

3.80 

Cl . 15) 

Gaussian 

A 

Mean 

-0.36 

C-0,11) 

-0,16 

C-0.049) 

-0.22 

C-0.06) 

Rayleigh 

*+ 

St, Deviation 

5,19 

Cl. 58) 

4.84 

Cl. 47) 

4.48 

Cl. 36) 

Model 

i ; 

Mean 

0.27 

CO, 08) 

-0,36 

C-o.ii) 

-0.20 

C-0.06) 

VLI 


St. Deviation 

3. 66 

Cl. 11) 

3.55 

Cl. 08) 

2.67 

CO . 81) 

Model 

c 

Mean 

0.20 

Co. 06) 

-0,10 

C-0.03) 

-0.33 

C-0.10) 

VLI 

o 

St, Deviation 

3.67 

Cl. 12) 

3.90 

Cl. 18) 

3.81 

Cl. 16) 

Model 
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TABLE 4 


NORMALIZED FOURTH AND SIXTH MOMENT DATA 
OF REAL AND SIMULATED TURBULENCE FIELDS 


Over a 10-min, sample) 


Model 

No. 


Real 
Atm . 

1 

2 

3 

4 

5 

6 


Normalized 

Gust 

Velocity 

Component 

Simulation 

Moment 

xa 

V 

w 

Technique 

Fourth 

3.5 

3.5 

3.5 

Real atmospheric 

Sixth 

21.7 

21.7 

21,7 

turbulence data 

Fourth 

3.0 

3,0 

3.0 


Sixth 

15.0 

15,0 

15.0 

Gaussian 

Fourth 

5.9 

3.5 

3.2 

Modified 

Sixth 

61.0 

22,3 

16,8 

Gaussian 

Fourth 

5.1 

3.2 

2.8 

Modified 

Sixth 

46.7 

18.9 

11,9 

Gaussian 

Fourth 

3.7 

3.2 

3.3 

Rayleigh 

Sixth 

21.7 

18.1 

19.9 

Model 

Fourth 

3.5 

3.2 

3.5 

VLI 

Sixth 

20.8 

16.0 

21,8 

Model 

Fourth 

3.1 

3.2 

3.9 

VLI 

Sixth 

14.0 

16.1 

21.5 

Model 





Probability % 




Probability % 
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Gust Velocity w ft/sec (m/sec) 
Figure 12 Probability Density of Simulated Field 








Model 6: VFI 

x Simulated Data 



3.75 (1.14) 7.75 (2.36) 11.75 (3.58) 

Gust Velocity w ft/sec (m/sec) 

Figure 15 Probability Density of Simulated Field 
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distribution is plotted on the same scale. It has been established (5) 
that real atmospheric turbulence exhibits a higher probability of both 
smaller and larger gust velocities compared to a Gaussian distribution. 
A careful study of the probability density of the simulated field 
reveals a higher probability of larger gust velocities compared to a 
Gaussian distribution, however the distributions, with the exception 
of Model 6, do not show higher probability of lower gust velocities. 
Power spectral densities of the simulated turbulence models are 
presented in Figures 16 to 21, The higher frequency components are 
compared with a line of slope -2 which is a characteristic of real 
atmospheric turbulence. The power spectrum in the entire frequency 
range within the limits of experimental error is in fairly good agree- 
ment with the assumed Dryden form (Equations 2,1,4 to 2,1.6), 

The patchiness of each of the models is plotted in Figures 22 to 
24, The derivative of vertical gust component is plotted illustrating 
a varying intensity of patchiness. Model 6 presents patchy character- 
istics which closely match real atmospheric turbulence. 

Element of surprise is tabulated in Table 5, At present there is 
no criterion available to either quantitatively measure this phenomenon 
or to establish a basis of comparison. In this report, "sudden jump" in 
the velocity field is used to describe element of surprise. 


Model 1: Gaussian 





Model 2: Modified Gaussian 
□ Longitudinal Component 
o Lateral Component 
X Vertical Component 





.Lateral 


Line of 
Slope -2 


Vertical 


-Longitudinal 



Frequency - CPS 


Figure 17 Power Spectral Density of Simulated Field 







Spectral Density (ft/sec) 2 /(rad/sec) 






Spectral Density (ft/sec) 2 /(rad/sec) 
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Mode 'I 6: VLI 



Figure 21 Power Spectral Density of Simulated Field 







TABLE 5 


FREQUENCY OF ELEMENT OF SURPRISE 
OF SIMULATED FIELD * 



Frequency of 

Element of 

Surprise 

Model 

u 

V 

v 

1 

0,03 

0,0 

0.0 

2 

0.07 

0.07 

0.0 

3 

0,03 

0,00 

0,0 

4 

0.0 

o.,o 

0,23 

5 

0,03 

0,03 

0.0 

6 

0.27 

0.40 

0.0 


*For a 3.5 ft/ sec jump in velocity field) 



SECTION V 


TEST PROGRAM 

This chapter describes the flight simulator experiment, including 
the details of the aircraft simulated, the flight simulator, and the 
pilot performance task. 

5.1 Simulated Aircraft 

The aircraft simulated is the Canadian deHavilland DHC-6 Twin 
Otter. This particular aircraft is chosen as representative of light- 
wing-loading S1QL aircraft. In addition, there are pilots available 
with flying experience in the Twin Otter who can validate the simula- 
tion. 

Aerodynamic and dynamic stability parameters are listed in 
References C5) and C8) . A summary is given in Table 6. 

5.2 Aircraft Simulator 1 

The Visual Motion Simulator (VMS) at the NASA Langley Research 
Center, a synergistic notion-base simulator with the basic interior and 
instrumentation of a jet transport cockpit (Figures 25), was 
employed in this study. A schematic diagram of the simulator, its con- 
trol system, and its data output capabilities is presented in Figure 
26 (8), A CDC-6600 digital computer, used exclusively to operate the 
real-time simulators, was programmed with the aircraft flight condi- 
tions, stability derivatives, six-degree-of-freedom differential equa- 
tions of motion, and a simulator washout routine. The program 

1 This description has been adopted from Reference (8), 
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TABLE 6 


AIRCRAFT PARAMETERS (REFERENCE 8) 


w 


u 0 = 


11500 lb (51152 N) 

256,67 ft/sec (78,2 m/sec) 
0.045 


I = 16900 slug-ft 2 (22907 kg m 2 } 

I = 27600 slug-ft 2 (37411 kg m 2 ) 

I = 40600 slug-ft 2 (55031 kg m 2 ) 

z 


°o = 


= -1,3 


420 ft 2 (39.0 m 2 ) 


C, = 


C„ = 


D 

Si 


a 

C 

m 

a 

C L. - 
a 


C 


c 


0.3818 

0,045 

0.035 

5.7295 

0.1432 

-1.9098 

1,52 


-0.89 

- 0.12 

0.1215 


c = 6,5 ft (1.98 m) 
b = 65 ft (19.8 in) 
Longitudinal Stability Coefficients 


C D. - 0 
a 

C = -5,9 

m. 

a 

C L = 5.504 

q 

Cd =0 

q 

C = -23.948 
m 

q 

C =-2C. 
x D 

u 

C = -2 C r 
z L 

u 


Lateral Derivatives 


C «* -0.1 
y P 

“ -0.5488 
P 

C = 0,006 
n 

P 


Control Derivatives 


C = C_ - C n 
x L D 
a 

C = -C T - 
2 L 

a a 

C = -C n 

x. D. 

a a 

C = -C 
z. L • 

a a 

C x ” _C D 

q q 

c = -c T 

z L 

q q 


= 0.5 


C z - 0.13 
r 


-0.1855 


n 


, i 


-- 0.39 
= 0.00348 
- 0,2055 


n P 


n, 


= 0.0398 


“ - 0,1 


m. 


= -1.79 


= 0.45 


0.01 


a o 



8 ) 


REPRO^IT x OF THE 

oewttaj, , is poor 
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Figure 26 Block Diagram of Motion-Base Simulator Apparatus (Ref. 8) 
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integrates the equations of motion 32 times a second. These values are 
used by the simulator washout routine to determine the position of the 
simulator legs and the dynamic characteristics of the hydraulic 
actuators. Since the simulator is not capable of producing the 
magnitude and the duration of displacements, velocities, and ac- 
celeration of the real aircraft, the washout routine appropriately 
scales down the predicted motions of the real airplane to values that 
the simulator can produce without exceeding any of its design limita- 
tions. The washout routine also attempts to drive the simulator legs 
back to their neutral position following a disturbance from equilibrium 
in anticipation of a future disturbance. A detailed description of the 
physical dimensions and the performance specifications of the VMS may 
be found in References (9) and (10) . 

The simulator is equipped to provide two sets of landing conditions. 
In the first set of conditions the pilot is not given any external or 
"out the window 1 ' cue. In the second set the pilot is shown a visual 
display which generates a realistic landing scene. The aircraft motion 
signals, through a feed back loop, run a video camera on a scale model’ 
of the airfield and its surroundings (Figure 27). In addition, the 
simulator is equipped with ILS (Instrument Landing System) Instrumentation 
which includes both a flight detector and raw data display. A 20° and 
30° flap configuration is provided on the simulator for the landing 
approach. With the motion-base and all of the features described 
above, the simulator provides a realistic representation of a Twin Otter 
in a landing configuration. 


A 



i-i i i i i J J 



Figure 27 Landing Scene 
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5 ,3 Pilot Task Performance 

Seven pilots experienced in civil, military and research flying 
were employed in the first part of the test program. Test runs, each 
of 8 to 10 minutes duration, for each of the six turbulence models were 
made in one pilot session. During separate sessions, some of the pilots 
repeated the models in random order. It was decided to have the pilots 
in a level flight, constant altitude tracking task with no visual or 
"out-of-window" cues in order not to introduce too many variables that 
might distract tne pilots from their primary objective of trying to 
distinguish differences between various turbulence models. After each 
run, the pilot was asked for his comment on turbulence through the use 
of a flight questionnaire (see Appendix B) . 

In the final part, the landing approach test program, eleven pilots 
were employed to fly the VMS. Three of the eleven pilots each had ex- 
perience of over 2,000 hours of flying in the Twin Otter. Test runs, 
each of 6 to 8 minutes duration, for each of the turbulence models, 
were made in one pilot session. During separate sessions, some of the 
pilots repeated the six models in random order. It was decided to have 
the pilots fly a constant altitude tracking task until the glide slope 
was intercepted. The pilots were then required to approach on a 6° 
glide slope using both the ILS and the visual display. The simulation 
is terminated at touch do™. After each run, the pilot, through a 
flight questionnaire (see Appendix B) was asked to estimate the turbulence 
intensity, realism, relative amplitude of aircraft motions in each of 
the six degrees of freedom, patchiness, workload, and to give a Cooper- 
Harper handling quality rating (Figure 28) (11) for the airplane 
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turbulence interaction. Additional, questions explored the bases for 
the pilot's judgments. In addition the pilots were also asked, to estimate 
the altitude, terrain, and atmospheric stability in relation to their 
flying experience. 

During each run, continuous strip-chart recordings were made dis- 
playing time histories of various aircraft parameters for later analysis. 
These include three linear accelerations , and three angular rates of 
the aircraft in body axes, elevator, aileron, and rudder deflections, 
throttle position, altitude, rate of climb and aircraft heading as well 
as the integral-squared error of the ILS track. A sample strip-chart 
is presented in Figures 29a and 29b. The nns intensities of the 
longitudinal, lateral and vertical gust fields are presented in Table 7. 

In the following chapter these opinion ratings will be statistically 
analyzed to establish the most realistic turbulence model and to identify 
the variables that critically affect the handling quality of aircraft 


in turbulence. 




HANDLING QUALITIES RATING SCALE 



ADEQUACY FOR SELECTED TASK OR 
REQUIRED OPERATION* 


AIRCRAFT 

CHARACTERISTICS 

DEMANDS ON THE PILOT 
IN SELECTED TASK OR REOUIREO OPERATION* 

PILOT 

RATING 

Excellent 
Highly desirable 

Pilot compensation not a lactor lor 
desired performance 

n 

Good 

Pilot compensation not a tacior for 


Negligible deficiencies 

desired performance 

Fair — Some mildly 

Minimal pilot compensation required lor 

El 

unpleasant deficiencies 

desired perlormance 



Is it 

satisfactory without 
. improvement? 


s' Is adequate^. 
s' performance 
^(tamable with a tolerable^ 
N. pilot workload? 


Deficiencies 

warrant 

improvement 


Deficiencies 

require 

improvement 


Minor but annoying 
deficiencies 

Moderately ob|ectionable 
deficiencies 

Very ob|ectionable but 
tolerable deficiencies 


Ma|or deficiencies 

Major deficiencies 
Majot deficiencies 


Desired performance requires moderate 
pilot compensation 

Adequale performance requires 
considerable pilot compensation 

Adequate performance requires extensive 
pilot compensation 

Adequate performance not attainable with 
maximum tolerable pilot compensation 
Controllability not in question 

Considerable pilot compensation is required 
for con'iol 

Intense pilot compensation is required to 
retain control 



§ 

o 

© H 
£5 B 
t=2 



it controllable ? 


Improvement 

mandatory 


Major deficiencies 


Control will be lost during some portion ol 
required operation 


Pilot decisions 


Cooper-Harper Rel NASA TNO-S153 


♦ Delimlion ol required apeiaiion involves oesigna'icn ol tliqhi phave and. or 
subphases wlh accompanying conomons 


Figure 28 Cooper-Harper Pilot Rating Scale 
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I 


a x longitudinal acceleration (g's) 

a y lateral acceleration (g's) 
a 2 vertical acceleration (g's) 
p roll rate (deg/sec) 

q pitch rate (deg/sec) 

Y yaw rate (deg/sec) 

blank 
time (sec) 



Figure 29a Typical Stlp Chart Output of Flight Parameters 


Ln 



3e Elevator deflection 

6a Aileron deflection 
6r Rudder deflection 

6t Throttle position {%) 
h Altitude (feet) 

Rate of climb (ft/mln) 
V Air speed (knots) 

T Heading (degrees) 


Figure 29b 
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Strip Chart Output of Flight Parameters 


TABLE 7 


INPUT-OUTPUT INTENSITIES 

Input rms ft/sec (m/sec) Output rms ft/sec (m/sec) 


iodel 
No , 

u 



V 



w 



ii 


v_ 


w 

1 

4,0 a. 3i) 


4,0 Cl, 31) 


4,5 (1,37) 

3,97 

(1.21) 

3,90 

(1.19) 

4,43 

(1.35) 

2 

m i= 3.1 

CO, 94) 

m 

- 3,2 

(0.97) 

m 

- 2,8 

(0,85) 

3,92 

(1.19) 

3.53 

(1.07) 

2,6 

(0,79) 


v 2 = 1.2 

CO, 37) 

V 

= 1.3 

(0,37) 

V 

- 0,9 

(0,27) 







3 

m - 3 , 2 

(0.97) 

m 

= 3,5 

(1.07) 

m 

- 4.1 

(1,25) 

3,90 

(1.19) 

3,90 

(1.19) 

3,80 

(1.16) 


v ” 0,8 

(0.24) 

V 

- 1,0 

(0,30) 

V 

= 0.9 

(0,27) 







4 

m - 6,1 

(1.86) 

m 

= 6,1 

(1,86) 

m 

= 6,1 

(1,86) 

5,19 

(1.58) 

4.84 

(1.47) 

4.48 

(1.36) 

5 

m = 3,1 

(0.94) 

m 

- 3,2 

(0,97) 

m 

= 2.8 

(0.85) 

3,60 

(1.09) 

3.5 

(1,07) 

2,60 

(0,79) 


v = 1. 2 

(0.37) 

V 

= 1.2 

(0,37) 

V 

= 0,9 

(0,27) 







6 

m = 3,2 

(0,97) 

m 

« 3.5 

(1.07) 

m 

= 4.1 

(1,25) 

3,60 

(1,09) 

3.92 

(1,19) 

3.82 

(1.16) 


v = 0.8 

(0,24) 

V 

« 1.0 

(0,30) 

V 

= 0.9 

(0,27) 








*m - mean 
2 v - variance 


SECTION VI 


RESULTS AND DISCUSSION OF SIMULATION 

Data obtained during the flight test program consisted of pilot 
opinion ratings and commentary ralating to the simulated environment, 
aircraft handling quality and data relating to the physical environment 
to which the pilot was exposed. The pilot opinion rating, obtained 
through a questionnaire, were in the following form: 

1. Realism of turbulence; 

2. Correctness of relative amplitude of disturbances; 

3. Patchy characteristics; 

4. Frequency content; 

5. Element of surprise; 

6. Atmospheric conditions; 

7. Handling quality ratings (Cooper-Harper); and 

8. Pilot task performance error. 

Since the flight experiment was conducted in two parts, the 
statistical analysis and other data will be presented separately. Figure 
and Table numbers with subscript A are the results of the constant altitude 
tracking task with no visual aid. In the same way subscript B denotes 
the results of the landing approach test program. 

The pilot opinion ratings have been statistically analyzed and the 
results are presented in the following forms: 

a) Mean and Standard Deviation : of pilot opinion ratings for each 

of the turbulences models (Figures 30A,B to 35A,B and 36A, 37A and 38B) . 

b) Correlation Matrix : Correlation among various physical 
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Figure 30B Pilot Opinion Ratings 
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Figure 31 A Pilot Opinion Ratings 
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Figure 32A Pilot Opinion Ratings 
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Figure 34A Pilot Opinion Ratings 
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Figure 34B Pilot Opinion Ratings 
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Figure 35B Pilot Opinion Ratings 
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Figure 36A Atmospheric Condition 
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characteristics of atmospheric turbulence is determined by 

a 

t (x. - x)(y - y) 

_ i=i 

Y xy " ri-1 CT x a y 

where 

Y correlation between x and y 

'xy 

x,y mean of x,y 

a ,a standard deviation of x,y 

x* y 

n number of observations. 

The correlation matrices are presented in Table 8A and 8B. 

The following observations can be made from the statistical 
analyses of pilot opinion ratings: 

1. Figures 30A and 30B present the pilot opinion ratings of handling 
quality and the realism of turbulence. It may be observed that at 
approximately the same rms intensity (see Table 7) of turbulence, the 
handling quality ratings transit from the satisfactory level, for a 
simple Gaussian model, to an unacceptable level for the more realistic 

and compositely structured VLI turbulence model. By comparing Figures 
30A and 30B it may be observed that visual display did not significantly 
alter the trend. 

2, Figures 31A and 3LB depict the element of surprise and the 
patchiness ratings. The Gaussian model (Model 1) was found to be a 
little too continu us by almost all the pilots under both experimental 


( 6 . 1 . 1 ) 

r 
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conditions. On the other hand, the Rayleigh model (Model 4) was rated 
"about right" as was the VLI model (Model 6) . 

3. Figures 32A and 32B present the frequency content (low and high) 
ratings. The Gaussian Model (Model 1) was poorly rated whereas the mean 
ratings of the Rayleigh (Model 4) and VLI (Models 5 and 6) turbulence 
models were in the range of "about right". In addition, the spread 
(standard deviation) in the pilot opinion ratings is greater during the 
visual landing approach task. 

4. Figures 32A,B to 35A,B present amplitude of disturbances as 
perceived by the pilots. The ratings show a progressive improvement as 
the pilots are exposed to more sophisticated models (see Models 4, 5 
and 6) . 

5. Figures 36A and 37A present the atmospheric condition observations 
in the form of terrain, altitude, and atmospheric stability. The primary 
purpose of evaluating these was to determine how sensitive the pilots 

were to changes in atmospheric condition. Most pilots, when exposed to 
the six turbulence models, thought they were flying over level plains. 

On the altitude rating, the pilots flying the Gaussian model felt this 
turbulence was typical of altitude greater than 10,000 feet whereas they 
consistently rated the other models as typical low altitude turbulence. 

6. Figure 38B presents pilot estimates of task performance. This 
is a measure of how well a pilot performed in tracking the ILS. Here 
we observe that the pilots had a greater difficulty in tracking the ILS 
for the VLI models than the simple Gaussian model. 

7. Table 8A and 8B present the correlation matrix for various 
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turbulence properties and aircraft handling qualities. Several important 
observations can be made from these symmetric matrices. From Table 8A 
(Level flight) , the realism of turbulence is highly correlated with 
patchiness (0.58), element of surprise (-0.63), and frequency content 
(0.52). This shows that in the opinion of pilots, the realism of a 
turbulence model is closely linked to the physical properties of real 
atmosphere. In addition, the high correlation between handling 
qualities and realism (0.74) indicates that the handling qualities are 
considerably worse for more realistic turbulence models. The low 
correlation between the patchiness character: sties and the intensity of 
turbulence (0.07) shows that the non-Gaussian patchiness characteristics 
cannot be induced by simply chosing a higher level of intensity (rms) . 

On the other hand patchiness is correlated to frequency content (0.45) 
and handling quality. Table 8B (the second part of the experiment) 
shows the same trend. Here pilot error is highly correlated with 
handling quality, and patchiness. 



TABLE 8A 


CORRELATION MATRIX 



Turbulence 

Intensity 

Realism 

Patchiness- 

Frequency 

Content 

Element of 
Surprise 

Handling 

Quality 

Turbulence Intensity 

1,0 

0,36 

- 0.07 

-0.16 

0,54 

0.27 

Realism 

0.36 

1,0 

0.58 

0.52 

-0.63 

0.74 

Patchiness 

0.07 

0,58 

1.0 

0.45 

0,39 

0.50 

Frequency Content 

-0.16 

0.52 

0,45 

1.0 

0.23 

0.63 

Element of Surprise 

0.54 

-0,63 

0,39 

0,23 

1,0 

-0,02 

Handling Quality 

0.27 

0.74 

0,50 

0.63 

-0.02 

1.0 



TABLE 8B 


CORRELATION MATRIX 



Turbulence 

Intensity 

Realism 

Patchiness 

Frequency 

Content 

Element of 
Surprise 

Handling 

Quality 

Pilot 

Error 

Turbulence Intensity 

1.0 

0.41 

0.21 

0.18 

0.27 

0.31 

.43 

Realism 

0.41 

1.0 

0.60 

0.49 

0.36 

0.70 

.38 

Patchiness 

0.21 

0.60 

1.0 

0.46 

0.38 

0.47 

.54 

Frequency Content 

0.18 

0.49 

0.46 

1.0 

0.31 

0.66 

0.11 

Element of Surprise 

0.27 

0.36 

0.38 

.31 

1.0 

0.31 

0.16 

Handling Quality 

0.31 

0.70 

0.47 

. 66 

0.31 

1.0 

.68 

Pilot Error 

0.43 

0.38 

0.54 

.11 

0.16 

0.68 

1.0 


03 

□ 



SECTION VII 


CONCLUSIONS 

This report has described several proposed turbulence models for 
producing artificial turbulence time histories which match the desired 
statistical properties of real atmosphere better than the presently- 
used simulation techniques. The use of these models gives improved 
realism and accuracy in piloted simulator studies of handling 
qualities as effected by atmospheric turbulence. 

From the analytical study of the time histories generated by these 
models, and their comparison with real atmospheric turbulence, the 
following conclusions can be drawn: 

a) Turbulence simulated by the VLI gust models adequately 
matches the probability distribution (fourth and sixth 
normalized moments, probability density, and cumulative 
probability) of real atmospheric turbulence; and hence, 
presents an improved representation of atmospheric 
turbulence, 

b) Frequency content and the patchy characteristics of real 
turbulence can be closely matched. 

c) The proposed turbulence models (VLI) can flexibly 
accommodate changes in atmospheric conditions characterized 
by terrain, altitude, and atmospheric stability. This 
flexibility is not provided by any of the presently-used 
techniques . 
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d) The mechanization of the proposed models on a motion-base 

simulator is easy and inexpensive computationally because these 
models utilize only three linear filters for the entire simulation. 

The time histories derived from turbulence models and the commonly 
used Gaussian model * were employed in a flight simulator experiment to 
determine the extent of pilot sensitivity to realism of various turbulence 
models and to evaluate the effects of turbulence on aircraft handling 
qualities. The principal conclusions drawn from the flight simulator 
study are: 

a) As expected from the analytical study, the pilot opinion ratings 
show a considerable improvement of i.u.^bulence properties (tealism, 
patchiness, frequency content, etc.) over the most commonly used Gaussian 
turbulence model. 

b) Aircraft handling quality and pilot task performance are critical- 
ly affected by low altitude clear air turbulence. Realistic turbulence 
models present greater difficulty in controlling the aircraft than 

simple Gaussian model. 

c) The correlation coefficient between the handling quality and 
the realism of turbulence is 0 . 74 . This high correlation indicates that 
the handling qualities are considerably -"urse for more realistic 
turbulence models. 

d) From the flight test results of this program, it is apparent 
that the pilot's ability to handle the airplane in a turbulent environ- 
ment not only depends on the rms intensity, but also the composition and 
the structure of turbulence. Pilots rated handling qualities in the 
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satisfactory range while flying in a turbulence environment simulated 
by a simple Gaussian model; whereas the handling quality ratings de- 
graded while flying in a turbulence environment simulated by the VLI 
turbulence model of approximately the same intensity. In fact, the 
handling quality ratings monotonically degrade as the pilots encountered 
more complex and realistic turbulence models. It may be concluded, 
therefore, that handling quality studies, using motion-base simulators, 
are critically affected by the suitable choice of a realistic turbulence 
model in addition to the appropriate rms intensities of turbulence. 

The tests were conducted in a simulated environment of a light 
general aviation STOL airplane. Caution should, therefore, be exercised 
in applying and extending the results to a general aircraft configuration. 
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APPENDIX A 


REVIEW OF BASIC DEFINITIONS C5) 

Stationarity : A random process is stationary if its statistical 

properties are not dependent on the time of their measurement. One 
could, for example, collect an infinite number of time histories, called 
an ensemble, which are representative of the process. If one takes an 
average across the ensemble, and if these averages are not a function 
of time, the process is stationary. 

Homogeneity : A random process is homogeneous if its statistical 

properties are independent of position, 

Ergodicity : In turbulence measurements it is impossible to obtain 

an ensemble from atmospheric measurements. Thus it is necessary to use 
time averages to get statistical information. If such a time average 
yields the same statistical properties as the ensemble average, the 
process is called ergodic. 

Mean Value : The mean value of a random variable, u, of an ergodic 

random process Is given by 


- = Lim _1_ rT 
u 2T J -T 




CA.l) 


In practice the limit Is not required and u can be approximated by 


u = f / 0 u(t)dt, for T large. 


(A. 2) 


This approximate representation is especially useful for processes such 


35 
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as turbulence. However, the time interval T must be large enough so 
that the average approaches the asymptotic value one would obtain for 
a stationary process. 

Variance : The variance of u is defined as 

^/! T CCuCt> -t)2]dt. CA.3) 

As before in practical applications the variance can be approximated 
by 

( . 

°u ^ T -^0 ” *hPdt, for sufficiently large T, (A, 4) 


Standard Deviation (Root Mean Square ): The standard deviation is 

defined as the square root of the variance. 

til 

Normalized Central Moment : The n normalized central moment, 

of a random process, u(t), is 


M 


n 


Lira i ,T - u d 

T-^° 2T J -T a 

u 


n = 1,2,3... 


(A. 5) 


•which can be approximated by 

M n " T ^0 [ ~ o n= 1,2,3... (A.6) 

u 

Cumulative Probability Distribution : The cumulative probability 

distribution of u(t) , F^Cx) is defined as the probability that u x. 
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Probability Density Distribution ! Probability density distribution 
of uCt), P^Cx) is defined as the probability that 

x < u x + dx, 

Gaussian Probability Density Distribution : If a random variable, 

uCt), is Gaussian distributed its probability density is given by 

F (x) exp C- i (~ — -) 2 D CA.7) 

U rz— 2 cr 

o /2ir u 

u 

Rayleigh Distribution : Another probability density of interest 

is the Rayleigh distribution defined as follows: 

2 

PCx) - — exp (- -1 — ) CA.8) 

c 2 2 c 2 

where c 2 is one half the expected value of the random variable x or 

c 2 « EGO = \ /q xP x (x)dx (A. 9) 

Cross Correlation Function : The cross correlation function of two 

random processes u(t) , wCt) is defined as 

R uw Ct) = jjf (t T uCOwCt + T)dt (A. 10) 

correlations are the measures of the predictability of a signal at some 
future time (t + x) based on the knowledge of a signal at time t. 


Autocorrelation Function: The autocorrelation function is a 


special case of the cross correlation function defined above in which 
wCt) = u(t), such that. 


R (t) 
uu 


Lim 


-ff f^ T utOutt + T)dt, 


(A. 11) 


Integral Scale Length : A statistical parameter of special 

importance in atmospheric turbulence is the integral scale length. 


L = — /“ R COdr, 

U n J -co UU 


U 


(A. 12) 


where is the reference steady state flight speed of the aircraft 
flying through turbulence. Scale length is an approximate measure of 
the distance an aircraft flies through turbulence. 

Cross Spectral Density : The cross spectral density of two random 

processes u(t) and w(t) is defined as the Fourier transform of their 
cross correlation 


* uw Cf> " /l R uw ^ T ^ ex P^" i2irf OdT , (A_. 13) 

where f is frequency. 

Rower Spectral Density : The power spectral density, PSD, of a 

random process is the Fourier transform of its autocorrelation function. 
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$ U C£) - /“ M R uu (T)expCi27rfT)dT. (A. 14) 

The PSD can be interpreted physically as the average contribution to 
the variable cr^from the frequency component f. Thus, 

a 2=3 J” <j> C£)df. (A. 15) 

U J -co '■ 

White Noise : White noise is a random process for which the PSD is 

a constant independent of frequency. That is, 

<f>^(f) - <|>q - constant. (A. 16) 



APPENDIX B 


FLIGHT QUESTIONNAIRE 

F] ight Number Date 

Pilot: 

1. Turbulence Intensity: 

'Light Moderate Severe Extreme 

2. Realism of Turbulence: 

Very Good Good Fair Poor Very Poor 

3. Correctness of Relative Amplitude of Disturbances: 

Not Enough About Right Too Much No Comments 

Ro I I 

Pitch 

Yaw ___ 

Heave 

Side Force 

4. Patchy Characteristics (Variation of Intensity Bursts) 

Much Too Continuous A Little Too Continuous About Right 

A Little Too Patchy No Comments 

5 . Frequency Contents of Turbulence: 

Not Enough About Right Too Much No Comments 

Low FRO: _ 

High FRQ: 


6. Element of Surprise 

In the Simulated 

Turbulence 

Field: 

a. Quite Often 

Sometimes 

Never 

— 

b. Real ism of 6a: 




Very Good 

Good Fair 

Poor 

Very Poor 
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7, Atmospheric Conditions: 

a. Altitude: 0 - 1,000 Ft 1,000 - 10,000 Ft 

Over 10,000 Ft Unable to Judge 

b. Atmospheric Stability: Stable Unstable 

Neutral Unable to Judge 

c. Terrain: Mountains Plains Unable to Judge 

8, Pilot Estimate of the Work Load: 

Very Easy Easy Average Difficult Very Difficult 

9, Pilot Estimate of Task Performance: (Integral Squared Error for ! LS 

Tracking Task) 

Very Good Good Average Poor Very Poor 

10. Realism of This Model Compared to Previously Flown Model: 

Very Good Good About the Same Poor Very Poor 

11. Did You Observe a Repetitive Pattern in the Turbulence Field?' 

Yes No 

12. Cooper-Harper Rating: 

13. Additional Comments About Realism of Turbulence and Aircraft Simulation: 
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APPENDIX B 
(Cont.) 

PILOT EXPERIENCE 

1 . Name Date 

2. What Type of Flying Experience Have You Had? 

Mi ! t fa ry __________ Civil 

3. Main Types of Aircrafts Flown: _ 


4. Total Number of Hours Flown: 

5. Hours of instrument Flying: ' 

G. Hours in Simulators: 

/. Hours in VMS: 

8. Hours in Twin Otter: 

9. a. Estimate the % of Time Flown in Turbulence: 
b. Of This Time What % Was Flown in 

Light Turbulence Moderate Turbulence Severe Turbulence Extreme Turbulence 


10. What Characteristic of Turbulence Interferes Most with Your Ability to 
Control the Aircraft? 


II. Describe the Most Critical Case of Turbulence Encountered During Your 
Flying Experience: 
a. Day Night 


b. 


Terra i n : 


Altitude*. 



Atmospheric Stability: 


Stable Neutral Unstable Unable to Judge 

What Was the Task You Were Attempting Before Turbulence Was Encountered 

(e.g. 1 LS Approach, Cruise, etc.) , 

Any Additional Comments: _ 


